Directional motion of forced polymer chains with hydrodynamic interaction 
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We study the propulsion of a one-dimensional (ID) polymer chain under sinusoidal external forces 
in the overdamped (low Reynolds number) regime. We show that, when hydrodynamical interactions 
are included, the polymer presents directional motion which depends on the phase differences of the 
external force applied along the chain. Moreover, the velocity shows a maximum as a function 
of the frequency. We discuss the relevance of all these results in light of recent nanotechnology 
experiments. 

PACS numbers: 36.20.-r, 47.15. G-, 62.25.-g, 87.19.ru, 87.10.-e 



I. INTRODUCTION 

Propulsion at the low Reynolds number regime is an in- 
teresting subject of investigation which induces us to re- 
view our intuition related to motion. The absence of iner- 
tial effects in this regime, along with the reduced number 
of degrees of freedom in small swimmers, makes the pos- 
sibility of movement not obvious at all. At this regime, 
the hydrodynamical interactions takes the essential im- 
portance being able to give the non reciprocity necessary 
to the net movement of a structured object. After the 
seminal work of Purcell [l| , and the one of Shapere and 
Wilczek [2] , an increasing interest has recently arisen and 
a good number of papers have appeared in the last decade 
to describe self propelled swimmers with, at least, two 
degrees of freedom, the minimal condition to guarantee 
the movement at low Reynolds numbers [![. The model 
proposed by Purcell, consisting of three joined rods, has 
been studied recently by Tarn et al. 3]. Other mod- 
els consider the non reciprocal movement of two spheres 
joined together by an extensible bond that are able to 
change their volumes Among the simplest models, 
one consisting of three spheres linearly joined together by 
means of extensible bonds [5|-[7| has been also proposed. 
Those one dimensional systems are actually the simplest 
possible devices with two degrees of freedom, where the 
motion is controlled by the phase difference between the 
relative motion of the spheres. An extension of the three 
(and more) particle model to two dimensions was studied 
by the authors of Ref . [|| . 

All these models deal mainly with the problem of self- 
propulsion, that is, the motion is caused byperiodic de- 
formation of their structure (see also Ref. Q for a com- 
prehensive review of these models), while in this paper 
we address a more general problem: the propulsion of a 
structured object (a polymer chain) submitted to non- 
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homogenous time dependent driving forces. Recently, 
experiments at nanometric scale have been performed by 
Dreyfus et al. @, Hj. In such works magnetic spheres 
joined by DNA links are driven by external fields. The 
velocity as a function of the frequency for this system 
shows a clear maximum. Moreover, a new experiment us- 
ing also oscillating magnetic fields driving a macroscopic 
extended body has been carried out by Garstecki et al. 
11], finding again a nonmonotonic dependence of the ve- 
locity versus the frequency of the applied field. Polymeric 
models that may qualitatively explain the experimental 
results have been implemented in Ref. [!, [l2j 

The system here studied presents rich and interesting 
properties. In fact, under this model, it is possible to ob- 
tain the above nonmonotonic behavior of the net velocity 
as a function of the frequency of the driving field. More- 
over, a strong dependence of the velocity with the phase 
difference of the force is present, showing that this phase 
represents a control parameter able to drive the motion 
in one direction or in the other, and also to modulate the 
maximum possible value of the velocity. 



II. MODEL 

Model equations. The model consists of N rigid cou- 
pled spheres joined to each other by a spring which can 
move in ID along its x axis (see Fig[T]). 

The system is supposed to be a few micrometers in 
size. By considering this scale and the typical velocities 
of microorganisms (or microscale or nanoscale devices), 
we obtain a Reynolds number around R ~ 10 -5 . For this 
reason, we can use overdamped dynamics in our equa- 
tions as a good approximation. At this scale we have 
to take into account the mutual interaction between par- 
ticles due to the fluid motion (i.e., the hydrodynamic 
interaction). For small particles moving in a fluid, the 
Oseen tensor approximation (13(, that relates the forces 
acting in each sphere with the respective velocities, can 
be used. It is worth noting that in overdamped dynamics 
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inertial effects are neglected and the system is not subject 
to acceleration. Consequently, in this approximation, the 
velocities in each instant depends on the instantaneous 
forces (through the Oseen tensor), and not on the past 
history of the system. 

Let Vi and F% be the velocity and force on the particle i, 
respectively. In one dimension, the Oseen tensor reduces 
to a square matrix with components ify which relates Vi 
to Ft, 
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The components of the Oseen matrix H are given by 
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where a and Xi stands for the radius and position of par- 
ticle i, respectively, and 77 is the fluid viscosity. 

In many of the previous works, and in particular in 
the three linked-spheres model [g, [7], the authors used 
internal forces to model the internal deformation of the 
system. Here, we drive the chain by an external force 
sinusoidal in time: 



Ki = Ai sin(wt - 



(3) 



terms represent, respectively, the elongation of the chain 
along the bonds of the sphere i with its neighbor. The 
parameter k is the elastic constant of the spring which 
describes the harmonic interaction between spheres. The 
hydrodynamical approximation under which the above 
equation maintains its validity is that at any time the 
distance between the spheres d is much higher than their 
radius a <C d. Thus, to fulfill this condition we used 
h = 40a > a. 

It is important to note that the external contribu- 
tion is in general different from zero [J^. Ki(t) 7^ 0], 
while the long time average of the external force is null 

[<Ei^i(*))t = o]._ 

The choice of a linear dependence instead of other more 
complex functions appears to be not crucial for obtain- 
ing the net movement given by our model. In fact we 
calculated the velocity with other distributions of phases 
and they gave the same qualitative dependence, and also 
the net flux for zero force as present in other studies in 
the literature [1]. We choose this linear distribution of 
the frequency essentially for two reasons: i) because it 
is simple and ii) because it is very easy to generalize it 
for an arbitrary number of spheres and to more complex 
(and realistic) systems of particles. 

For calculation purposes, we use an adimensional set 
of equations, where the length unit is the radius of one 
sphere (a), the time unit is 6irr]a/k, and the forces are 
measured in units of ka. The equations of motion ob- 
tained are then: 



where the amplitudes Ai are constant and equal for all 
the particles Ai = A, with i = 0, (N — 1) and <f> is the 
total phase difference between the first and the last par- 
ticle. The intermediate phases between the spheres are 
equally distributed along the chain, so that <f>i — i 
with i = 0, (N - 1). 

Therefore, our system has three free parameters: the 
amplitude A, the frequency u of the driving forces and 
the total phase difference (j> along the chain. 

The total force acting on the i-th particle, if we exclude 
the hydrodynamical contributions, is then: 

Fi = Ai8m(ut+4>i)+kAl i+ -kAli- i G [0, ...,N-1] (4) 



N-l 

Xi{t) = Vi = Kij(t) F i(t)> 
3=0 



where the new mobility matrix is given by: 



where we define Al i+ as 



l and AZj_ as x 



X j 



Iq, Iq being the natural length of the link. Both 
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FIG. 1: (Color online) Scheme of the chain compound 
by N monomers. The forces act on every sphere and 
the last one is shifted in phase with respect to the first 
one by a phase <p. The i-th sphere presents a phase shift 
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(6) 



The system is not solvable analytically for an arbitrary 
number N of particles. We solve numerically the equa- 
tions using a Runge-Kutta algorithm of the fourth order, 
by using a variable increment dt, adjusted for the differ- 
ent frequencies of the driving force. 

The main measure under study here is the asymptotic 
velocity of the polymer center of mass, and more precisely 
the mean velocity over a long time 7m, which is much 
bigger than the period of the applied force T = 2-k/uj. 



v = lim 

T M ^ NT M J t 



t +T M N ~ 1 
3=0 



(7) 



where to is here a proper transient time in the dynam- 
ics. This velocity depends on the total phase difference 
4> and on the frequency u> of the applied force. Although 
we have not been able to demonstrate the unicity of the 
velocity, the simulations (see below) indicate that this is 
the case independently on the initial conditions. 
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Model symmetries. The system under study presents 
interesting symmetries useful to demonstrate some of the 
features found in the numerical simulations. We first note 
that it is simple to see that the equations of motion are 
periodic in the total phase difference with period 4>p — 
(A — l)27r. In fact, every sphere is submitted to a force 
Asm(ujt + 4>i), and so, if we add a phase <pp to the total 
phase difference, (i.e., — > + 4>p), we obtain 

Asm(ujt + <f>i) — ► Asin(u>t + ijr^r + 27ri) 
= Asm{ujt + 4>i + 2-ki) = Asm(ujt + 

valid for any <f>i. The first property is then: 

v{^ + <f> P )=v{4>). (8) 

On the other hand, we can also see that with the trans- 
formation 4>i — > — <f>i the velocity changes sign. In fact, 
by using the time transformation the t — > — t ± T/2, we 
have 



A sin(wi + (pi) - 
= —A s'm(ujt - 



Asm(-Lut±uT/2 + = 
<pi ± 7r) = Asin(wi — <pi). 



This transformation maps the equation of motion for par- 
ticles with phase <fi to that one with phase —<j) and op- 
posite velocity. In fact, Eq. ([5]) can be formally written 
as: 



^ = f[xi,sm(ut + fc)]. 



(9) 



Using the above transformation in time (t — >• — t ± T/2), 
and the demonstrated sine transformation, the trans- 
formed equation reads: 



^ = f[xi,sin(ut- fa)]. 



(10) 



This equation is identical to the original one except 
for a global sign and for a change of the sign in all the 
phases. In other words, the second property is: 



v{4>) = v{-4>). 



(11) 



Note that this property allows us to choose the direc- 
tion of the velocity by changing the relative phase be- 
tween particles of the external driving forces. 

By using the property it is immediate to state 
that v(0) = 0. Remembering also the property ©, it is 
also immediate that v(4>p) = 0. Moreover, by using the 
property ([8]) with <j) = —<j)p/2, we have: 

= v(-4> P /2). (12) 

But using again the property (TTT|) it is also true that 

v{<j> P /2) = -v{-4>p/2). (13) 

So it comes out that v(c/>p/2) = 0. 

In synthesis, we can summarize these last outcomes as 
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FIG. 2: (Color online) Frequency (w) and phase (</)) 
dependence of the velocity for a different number of 
particles N. Note that for reading purposes, the panels 
with N — 5 and 10, have the scale of frequencies 
different from the others. 



III. 



SIMULATIONS. 



v(0) = v{<t) P /2) = v((j> P ) = 0. 



(14) 



We have performed extensive simulations with differ- 
ent numbers of spheres of the polymer. The main mea- 
sure is the asymptotic velocity (more precisely, the mean 
velocity over many periods of time) for different values 
of the driving frequency (uS) and for different values of 
the total phase difference between the last and the first 
particle (0). Amplitude is fixed to A — 10 along all the 
simulations. In Fig. [5] it is summarized the behavior of 
the asymptotic time average velocity as a function of the 
two variables say frequency u> and phase difference (j>, for 
different values of the chain length, namely N = 2,3,5, 
and 10. We can observe clear maxima and minima of the 
velocity which strongly depend on the number of particles 
N . In particular, by increasing A, the oj-<p region where 
the velocity is significantly different from zero reduces 
rapidly with the frequency, while the absolute values of 
the velocity maxima increases. These features are better 
illustrated in Figs. [3] and 03 where a few cross sections of 
the previous plots have been depicted. In Fig. [3] we can 
see that for N = 2 ((a) panel) a minimum and a maxi- 
mum appear for the velocities as a function of the phase 
difference, for fixed values of the frequency. In this case 
the periodicity in the total phase difference is given by 
4>p(N = 2) = 2ir. The case A = 3 seems more complex, 
but still with one minimum and one maximum, with a 
total phase difference 4>p(N = 3) = 4ir. In general the 
period depends on A: cj>p = (A — l)27r as it has been 
shown above. For clarity, panels (c) and (d) only show a 
semi-period. The three frequencies used in all the plots 
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(C) TV: 



(d) N = 10 



FIG. 3: (Color online) Phase difference ((f)) dependence 

of the velocity. The vertical lines represent the first 
phase difference at which the external total force -Ftot is 
at any time. Note that for N = 5 and 10 only a half 
period in phase is shown. 



(see the legend in each panel), have been chosen before, 
after and around the maxima values visible in Fig. [5J 

The vertical lines present in the Fig. [3] indicates the 
first phase values at which the total external force on 
the system is null at any time t. In that case the move- 
ment can be associated with an internal movement only, 
typically used for swimming, where the external force is 
zero. It is there visible that a net velocity is found for 
those values for N ^ 3, as also found in other works on 
swimming systems. Applying the sum rule for the sine 
function, the null external force is recovered for the phase 
value (f>= 2n{N - 1)/N. 

To better visualize the kinetic behavior of the chain, a 
set of trajectories of the center of mass have been drawn 
in Fig. @] for the case N = 3. We can observe in the 
upper panel that for <f> — 2n no displacement is present. 
The cases <j) = 2.0 and <j) = 10.28 represent respectively 
a value around the minimum and the maximum of the 
velocity, as visible in the corresponding plot of Fig. [3] 
The other phase value <j> — 4.29 lies close to the null ex- 
ternal force. In that case we can see that the center of 
mass does not present visible oscillations (because the 
oscillating external forces compensate) and the net mo- 
tion has a negative velocity, as also shown in Fig. [3] with 
N = 3 at the velocity given at the vertical line. The 
lower panel of Fig. @] shows the motion of each of the 
three spheres for the phase value <\> — 2.0. We can ob- 
serve the broken symmetry in the amplitude of the three 
different oscillations. The leading particle undergoes the 
larger oscillations. This is true for all values of the total 
phase difference independently of the velocity direction. 



80 
70 
60 
J50 
40 
30 
20 



_ i, - 2.00 



0.4.19 



600 

80 
60 
40 
20 


-20 g 
600 



650 



700 
t 



750 



800 



WVWWVWWWiWAMAAA/WsAM 



4> = 2.00 



650 



700 
t 



750 



800 



FIG. 4: (Color online) Upper panel: Trajectories of the 
center of mass of the polymer for the case N = 3 and 
for different values of the total phase <j). Lower panel: 
Trajectory of each of the three particles for the phase 
value cf> = 2.0. All trajectories are calculated with the 
frequency U) = 2.0. 



Figure [5] shows the cross section of the velocity as a 
function on the frequency, for different values of the phase 
difference. The interesting result is that a clear maxi- 
mum appears for all the chain lengths. This maximum 
decreases in the frequency value, as we increase the num- 
ber of spheres in the chain. In the extreme case of a very 
large number of spheres, we can argue that no movement 
persists, except at uj = (i.e. for a force constant in 
time). 

It is interesting to note that the null value of the veloc- 
ity corresponds to the maximum force amplitude to the 
system. In fact, for example, for <fi = 0, the force is the 
same for all the spheres, and the total external force on 
the chain is simply given by F to t = NA, where A is the 
amplitude of the force on each sphere. The reason why 
the mean velocity is zero, is that the system is in that 
case strongly symmetric, and the average over many time 
averages gives a zero mean velocity. 

Figure [S] shows the behavior of the values of the max- 
ima of the velocity, where we observe their rapid increas- 
ing trend as a function of the number N of spheres in the 
chain. This is due to the fact that the maximum veloc- 
ity occurs for a net force different from zero, and so an 
increase of the number of particles results in an increase 
of the value of the force acting on the chain. The insets 
of Fig. [5] show the frequency w max and the total phase 
difference </> ma x corresponding to the maximum velocity 
as a function of the number of spheres N. By increasing 
N, CJ max moves toward the zero value and a power-law 
scaling cj max s» N~ ls is also observed. 

Finally, we should stress that this non-monotonic be- 
havior of v with uj does not appear in the study of self- 
propulsion system like those of Refs. [H,0] where a simple 
v oc ui is observed. However, it is recovered in Ref. Q 
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FIG. 5: (Color online) Frequency (w) dependence of the 
velocity for different number of particles and for three 
different phases. 



when some driving force is applied. This behavior has 
been observed in the two experiments pil [H| cited be- 
fore on devices (at very different length scales) under 
external fields. 
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FIG. 6: (Color online) Modulus of the chain velocity 
maximum as a function of the number of spheres N. 
Insets: Frequency of the maximum in linear plot (top 
left) and log-log plot (top right) , and phase of the 
maximum (bottom). 



IV. DISCUSSION AND CONCLUDING 
REMARKS 

This work presents a model of directional motion based 
on the effect of modulated driving forces and hydrody- 
namical interaction between particles. The introduction 
of inversion-symmetry-breaking fields to provoke unidi- 
rectional motion has been suggested previously in other 
contexts, like the motion of molecular motors [lil . 
motion of solitons (la |. or active Brownian motion 
|l9( . In these cases it is needed either to have a substrate 
potential or inertial effects or the presence of fluctuations. 
On the contrary, this model deals with a chain in an over- 
damped environment without any external potential. 

It is interesting to stress that the net velocity is given 
by the hydrodynamical contribution, as it is simple to 
verify. In fact, if the friction is local, (i.e., the Oseen 
matrix of Eq. @ remains with the diagonal terms only) 
the equation of motion ([5]) becomes linear, and the sum 
of the forces (internal and external) over all particles, 
averages exactly to zero. In this sense, it is clear that 
the hydrodynamical interactions introduce a non local 
and non linear contribution which generates the peculiar 
net velocity here found. Therefore, this contribution is an 
essential term to be considered in these kinds of problems. 

Another interesting point to discuss is the fact that, in 
our system, the two particle chain presents net motion for 
a large range of parameters. This could appear in contra- 
diction with the so called 'scallop theorem" [HH3], which 
requires a system with, at least, two degrees of freedom 
to support the motion. However, one has to note that our 
system does not fulfill the other requirement of the 'scal- 
lop theorem". In fact, in that case, deformation of the 
swimming object should be driven by internal forces that 
average to zero at any time. This only occurs here for cer- 
tain total phase differences (f>, whose values are marked 
with a vertical line in Fig. EI In fact; for TV = 2, the theo- 
rem hypothesis (null force) is fulfilled for <f> = vr, and the 
velocity is zero at this point. Larger polymers (N > 2) 
show net velocity at these special null-force points. 

In summary, we analyze a polymer chain which moves 
in a low Reynolds number regime under a driving force 
which depends on the positions of the particles it is com- 
posed of. The study is mainly devoted to the influence 
of the driving force to the possible presence of the net 
velocity in the system. The presence of some symmetries 
in the equations has been analyzed as well as the rela- 
tionship of the chain velocity with the total phase differ- 
ence between the last and the first particles of the chain, 
finding a maximum, and revealing that the directionality 
can be changed by controlling the phase differences of 
the driving force across the chain. The maximum of the 
velocity is also present as a function of the frequency of 
the applied force, showing again an optimal value for the 
transport which depends on the number of particles in 
the chain. The model here presented includes a special 
case of self propelled swimming, given specifically by a 
total phase difference such that the total external force 
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acting on the chain be null. In that case the results are 
in agreement with other studies. The velocity behavior 
predicted by this model is also in agreement with the ex- 
perimental outcomes of time dependent driving forces on 
extended bodies. 

This study is a first approach to the effect of both a 
time and a spatial dependent field acting on elastically 
deformable objects. The driving proposed here can be 
performed in devices at microscopic scale using electro- 
magnetic waves where phase differences are related to 
the particle position. This means that it could be pos- 
sible to use a suitable light beam to drive the polymer. 
Our model is a simple approximation to this experimen- 
tal proposal if the particle oscillations are small in com- 
parison to the distance among them. The possibility of 



tuning the value and direction of the velocity with few 
parameters makes this system feasible to pilot structures 
in viscous environments. 
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